Graduated optimization

Graduated optimization is a global optimization technique that attempts to solve a difficult optimization problem by initially solving a greatly simplified problem, and progressively transforming that problem (while optimizing) until it is equivalent to the difficult optimization problem.[1][2]

Contents

Technique description

Graduated optimization is an improvement to hill climbing that enables a hill climber to avoid settling into local optima. It breaks a difficult optimization problem into a sequence of optimization problems, such that the first problem in the sequence is convex (or nearly convex), the solution to each problem gives a good starting point to the next problem in the sequence, and the last problem in the sequence is the difficult optimization problem that it ultimately seeks to solve. Often, graduated optimization gives better results than simple hill climbing. Further, when certain conditions exist, it can guarantees to find an optimal solution to the final problem in the sequence. These conditions are:

It can be shown inductively that if these conditions are met, then a hill climber will arrive at the global optimum for the difficult problem. Unfortunately, it can be difficult to find a sequence of optimization problems that meet these conditions. Often, graduated optimization yields good results even when the sequence of problems cannot be proven to strictly meet all of these conditions.

Some examples

Graduated optimization is commonly used in image processing for locating objects within a larger image. This problem can be made to be more convex by blurring the images. Thus, objects can be found by first searching the most-blurred image, then starting at that point and searching within a less-blurred image, and continuing in this manner until the object is located with precision in the original sharp image. (Representing images with varying levels of blurring is called pyramid representation, but pyramid representation is also used for other purposes besides finding objects with graduated optimization.[3])

Graduated optimization can be used in manifold learning. The Manifold Sculpting algorithm, for example, uses graduated optimization to seek a manifold embedding for non-linear dimensionality reduction.[4] It gradually scales variance out of extra dimensions within a data set while optimizing in the remaining dimensions. It, therefore, begins in a globally optimal state with a trivial problem, and as it scales variance out of the extra dimensions, it gradually transforms the problem into that of projecting data into a small number of dimensions.

It has also been used to calculate conditions for fractionation with tumors,[5] for object tracking in computer vision,[6] and other purposes.

Related optimization techniques

Simulated annealing is closely related to graduated optimization. Instead of smoothing the function over which it is optimizing, simulated annealing randomly perturbs the current solution by a decaying amount, which may have a similar effect.

References

  1. ^ Thacker, Neil; Cootes, Tim (1996). "Graduated Non-Convexity and Multi-Resolution Optimization Methods". Vision Through Optimization. http://homepages.inf.ed.ac.uk/rbf/CVonline/LOCAL_COPIES/BMVA96Tut/node29.html. 
  2. ^ Blake, Andrew; Zisserman, Andrew (1987). Visual Reconstruction. MIT Press. ISBN 0-262-02271-0. http://research.microsoft.com/en-us/um/people/ablake/papers/VisualReconstruction/. 
  3. ^ Crowley, James (1981). A representation for visual information (PhD thesis). Pittsburgh, Pennsylvania: Carnegie-Mellon University, Robotics Institute. OCLC 227552686. http://www.ri.cmu.edu/pub_files/pub3/crowley_james_l_1981_1/crowley_james_l_1981_1.pdf. 
  4. ^ Gashler, M.; Ventura, D.; Martinez, T. (2008). "Iterative Non-linear Dimensionality Reduction with Manifold Sculpting". In Platt, J. C.; Koller, D.; Singer, Y. et al.. Advances in Neural Information Processing Systems 20. Cambridge, MA: MIT Press. pp. 513–20. http://axon.cs.byu.edu/papers/gashler2007nips.pdf. 
  5. ^ Afanasjev, BP; Akimov, AA; Kozlov, AP; Berkovic, AN (1989). "Graduated optimization of fractionation using a 2-component model". Radiobiologia, radiotherapia 30 (2): 131–5. PMID 2748803. 
  6. ^ Ming Ye; Haralick, R.M.; Shapiro, L.G. (2003). "Estimating piecewise-smooth optical flow with global matching and graduated optimization". IEEE Transactions on Pattern Analysis and Machine Intelligence 25 (12): 1625–30. doi:10.1109/TPAMI.2003.1251156.